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Abstract 

A new model of physics, with a hidden conformal sector which manifests itself as an unparticle 
coupling to Standard Model particles effectively through higher dimensional operators, predicts 
strong collider signals due to unparticle self- interactions. We perform a complete analysis of the 
most spectacular of these signals at the hadron collider, pp{p) 7777 and 'j^gg. These processes 
can go through the three-point unparticle self interactions as well as through some s and t channel 
diagrams with one and/or two unparticle exchanges. We study the contributions of individual 
diagrams classified with respect to the number of unparticle exchanges and discuss their effect 
on the cross sections at the Tevatron and the LHC. We also restrict the Tevatron bound on the 
unknown coefficient of the three-point unparticle correlator. With the availability of data from 
Tevatron, and the advent of the data emerging from the LHC, these interactions can provide a 
clear and strong indication of unparticle physics and distinguish this model from other beyond the 
standard model scenarios. 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) at CERN, once functional, will attempt to discover the 
last piece of the Standard Model (SM) puzzle (the Higgs boson, responsible for electro-weak 
breaking) as well as hopefully provide signals of physics beyond the standard model (BSM). 
The main process responsible for producing a Higgs boson at the LHC is gluon fusion, and 
a light boson produced in such fashion is then expected to decay into two photons. Such 
a decay is supposed to be weak, as Higgs interact only weakly and in the SM the process 
occurs at one-loop level only. 

In a BSM scenario, a new scalar field which is singlet under the SM gauge group can 
mix with the Higgs boson but also couple to photons and gluons directly through higher 
dimensional operators with a cutoff scale A.u- Such a model has been proposed by Georgi l|, 
based on the hypothesis that there could be an exact scale-invariant hidden sector at a high 
energy scale. Below A^^, the model emerges as interpolating fields with general non-integral 
scaling dimension dy-, which behave like a du number of invisible massless particles-hence 
the term unparticle iU) used by Georgi. An unparticle does not posses a fixed invariant 
mass squared, but has a continuum mass spectrum, as a consequence of scale invariance 

p(/) = A,, 9{p') 9{p') (1) 

with Arf^ the normalization factor. The production of unparticles at low energies is described 
by an effective field theory and can give rise to peculiar energy distributions because of the 
non- integral values of du- 

While the concept dates only to 2007 there have been over 120 studies so far for 

n n 

implications of unparticle physics topics as varied as on neutrino physics [2[ cosmology [3|, 
rare processes in the standard model 4|, effects on precision measurements j^, as well as 
specific low energy signals for the detection of unparticles jo]. The peculiarities of unparticle 
physics have motivated numerous studies of their implications at the colliders 0, S]. Re- 
cently, it has been pointed out j^, [l^ that among collider signals, unparticle self-interactions 
give rise to spectacular effects in gluon fusion processes. (Note also the more recent study 
on unparticle self interactions [ill]). In particular, while the interactions gg ^ U ^ 'yy 
lead to enhanced signals in the Higgs boson decay channels, the three-point self-interactions 
gg ^ U UU 7777 are practically background free and may provide an unsuppressed 
and extremely promising signal for unparticle discovery. The unusual feature here is that 



the production of additional unparticles with high px in the final states does not suppress 



the production rate, unlike the production of known particles 



9| . There is encouraging data 



from the L3 collaboration on 77 interaction, which exceeds the prediction of the standard 



model by about one order of magnitude at the highest transverse momentum data [12|, and 
gluon fusion might show similar enhancements. Motivated by these considerations, in the 
present paper we extend the calculation in 9| and we present a complete evaluation of the 
relevant diagrams, in a comparative fashion, for the processes ppij)) 7777 and 77(7(7 at 
both the Tevatron and the LHC. We restrict our analysis to the hadron colliders: the LHC 
because of the large center-of-mass energy, and the Tevatron, because it operates now and 
could produce visible signals of physics BSM. A similar analysis is possible for ILC, but for 
our purposes, it is beyond the scope of this work. We also update the Tevatron bound on 
the three-point unparticle self interaction discussed in , based on the complete evaluation 
of the cross section. 

We proceed as follows: in Section Ull we describe briefiy the unparticle model features, 
concentrating on the ones relevant for our considerations. In Section IIIII we present the 
numerical investigation of processes generated by unparticle self interactions in 7777 and 
77(7(7 at both the Tevatron and the LHC. We compare the results obtained from production 
subprocesses, grouped by number of unparticles in intermediate states, and by channels. 
There are basically only two important unknown parameters, namely the scaling dimension 
du and the cutoff scale An. We restrict a third parameter, emerging from the three-point 
unparticle correlation function from the Tevatron data. We discuss the results obtained and 
conclude in Section IIVI In the Appendix we give the complete list of relevant diagrams for 
our calculation. 



II. UNPARTICLE FORMULATION 

In this section, we summarize the basic features of unparticle physics. In order to be 
relevant, unparticles must interact with SM fields, although the detailed dynamics of the 
interactions are not known. These interactions are parameterized in the effective Lagrangian 
approach as interactions between a new physics operator Ou with dimension du and the 
standard model operator Osm with dimension dsM 
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with \dsM ^ coupling constant of order 1. From scale invariance considerations, the phase 
space for the unparticle operator is a continuous spectral density 

dusF,, = A,, e{p') e{p') {pr^-'-0^, (3) 

where is the normalization factor 

A - ^^^^ + 

i2nr^u r {du - 1) r i2du) ^ ^ 

Unparticle operators interact with the SM fields through the exchange of heavy particles of 
mass M. Integrating out the heavy fields produces a series of effective operators describing 
unparticle interactions with the SM fields at__low energy The operators describing the 



interactions for various unparticle spins are 



SpinO: Kj^JfO, K-^J^^'fO, X.^G^.G^^O 
Spin 1 : Ai-^ 77,/ , Al-i^T 77,75/ , ■ ■ ■ 
Spin 2 : -^K^i^^ (l, +7. D,) , X^-^G.^G^^O'^ 

For the spin-0 case, the interaction of U with a fermion pair is suppressed if there is a 75 in 
the vertex and we neglect such terms. Note also that, unlike the above form, in our study 
this vertex is defined with a Higgs coupling such that the prefactor becomes eX'^v/ A.'jj where 
V is the Higgs vacuum expectation value. Otherwise, one needs to introduce two cutoff scale. 
The propagators for the unparticles were suggested to be j^: 

^^^^'^ = 2sintrf«(|p2|+Io+)2-'^« "^'^ 

for the scalar, and 
and 

for the spin-1 and spin-2 unparticle, respectively. The phase is defined as = Arg [(— p^)"^"] . 
Thus, it is non-zero only in the s-channel where is positive while is identical to zero in 
the t— and u— channels. Of course, the +zO"'" piece is only relevant in the s channel. Here 



For a scalar operator, unitarity requires as a lower bound, (4/ > 1. While no upper 



bound exists, supersymmetric examples suggest values of du < 2 



13| . For vector unparticle 



operators du > 3 and for symmetric and antisymmetric tensor operators du > 2 and du > 4 
respectively [l^. The cross section is dependent on du as and it becomes much 

smaller for the values of du larger than three. For this reason in this work we only consider 
the scalar unparticle effects. 

A scalar unparticle couples to the Standard Model sector through: 



■u 



(8) 



The Feynman rules for the scalar and tensor unparticle operators coupling to the gg and 77 
are, respectively 

U 



O : 



9,1 ^du 

u. 



(9) 



where K^^f"' = —g^'^PiP'^ ~ 9^^9^'^Pi ■ P2 + 9''^PiP2i A°, are the scalar and tensor 
coupling constants. 

Evaluation of the processes ligg, 7777 require evaluation of both the two point 



correlation function 



{0\Ou{x)Ol{Q)\0) 



d^P 



-iPx 



as well as the three point interaction 



9| 



\Ou{x)Ou{y)Ol{m 



(|x - y\ \x\ \y\Y^ 

The later yields, in momentum space, using the two point correlation function: 

{Q\Ou{Pl)Ou{p2)OUPl+P2m 



(10) 



(11) 



dS 



{ [-q^ - «e] [-{Pi - q? - ^e] [-{P2 + qf - le] ] 



n-2 



Pi P| 

s s 



(12) 



where n = Q{1 — du/4:) and 
T{n-2) 



9 2 
Pi P2 



16vr2 [rm'Jo 



dxidx2dxs{xiX2X3) 3 ^6{xi + X2 + X3 — 1) 



n~2 



(13) 



with A = xiXovi/s + xix^p\/s + and s = {pi + ^2)^- As a check, we reproduced 



the Fig. 2 in 
agreement. 



9| where Fy is plotted as a function of its one argument and got complete 



In all the analyses that follow we take A^^ = 1 and Ag = v27r/e. We restrict our 
considerations to the scalar unparticle, as it was shown that the vector and tensor give 
smaller contributions 



III. NUMERICAL ANALYSIS 

The complete set of Feynman diagrams for the subprocesses gg, qq ^ U... 7777 and 
for the subprocesses gg, qq U... 11 gg are given in the Appendix. There are quite 
a number of diagrams contributing to the processes, requiring some automatization to do 
calculation. We generated our results using the softwares FeynArts and FormCalc jl^ and 
the Vegas Monte Carlo program for the four-body phase space integration as well as for 
the numerical evaluation of the function Fy. In order to avoid coUinear singularities, we 
employ angular cuts as well as cuts on the transverse energies of the final photons and glu- 
ons. We require the final photons and gluons to have at least 30 GeV transverse energy 
and their ejection directions make at least 10 degrees with the beamline in the forward and 
backward directions. Since the unparticle model has unconventional features, the vertices 
and propagators the softwares FeynArts and FormCalc required some modifications. The 



hadronization is done by using the LHAPDF version 5.3 [17| with the parton distribution 
functions of CTEQ6LL. We neglect the masses of the light quarks in the proton while calcu- 
lating the contributions from the quark pair. We also neglect the contributions from the 
SM diagrams but include the SM background separately. In fact, the SM backgrounds for 
PP — > 7777 and the signals pp 1199 go only through quark subprocesses at partonic level. 
Under the same conditions, we estimated the background cross sections as follows: at the 
Tevatron; 2.6 x 10~^ pb and 0.3 pb for pp — > mi and pp ligg , respectively; at the 
LHC, 7.7 X 10"^ pb (for pp — > mi ) and 0.8 pb (for pp ligg )• 

In order to keep track of various contributions and make comparisons, we group the 
graphs according to the partonic constituents {gg or qq), the number of unparticles present 



Note that the vector unparticle does couples only qq while the tensor couple to both qq and gg. 
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in the intermediate states, and whether they proceed through the s or the t channeL We 
follow the same path in the numerical evaluation of each contribution, which we discuss in 
the following subsections. 

The existence of more than one external gluon requires an extra care. To preserve the 
gauge invariance one either needs to include the ghost contributions or perform the the 
gluon polarization sum for only the transversal polarization. We follow the second path by 
introducing appropriate projection operators for the transverse polarization states of the 
gluon. Since one can still do the sum conventionally for one of the gluons, in the gg — > 7777 
and gg ligg subprocesses, we need a single and three projection operators, respectively, 
when summation over the gluon polarizations is performed (for more details, see |l8l|). 



A. The multiphoton processes at Tevatron 

While the LHC seems to be still besieged by technical problems, it is worthwhile to 
explore signals coming the physics of the currently running largest energy collider facility, 
the Tevatron at Fermilab. Tevatron has accumulated yet unexplored data, and can still 
expected to yield signals for physics beyond the standard model. With this in mind, we 
begin the exploration of unparticle interactions in the multiphoton signals at the Tevatron, 
specifically looking at the signals pp — > 7777 and pp ^^gg ■ The center of mass energy 
^/s = 1.96 GeV for Tevatron is used throughout. 



1. The process pp 7777 at Tevatron 

We explore first the dependence of the total cross section on the value of the parameter 
A.U- Fig. [1] shows the variation of the cross section with the unparticle energy scale for three 
values of dy (1.1, 1.5 and 1.9). As expected, the cross section is the largest for the smallest 
du- Within each panel, we plot separately contributions from the two partonic channels 
{gg and qq) and momentum transfer channels {s and t), and present separate contributions 
coming from one, two and three unparticle exchange. One of the processes, the s-channel 



contribution coming 
qq was calculated in 



Tom the exchange of three unparticles as intermediate states in gg and 
9|. Clearly, as can be seen from all three graphs in Fig. [H this process 



is subdominant for all values of A^. As expected, the processes with the least number of 
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unparticles dominate, as we must include in the evaluation one order of the energy scale Ku 
for every unparticle exchanged. In addition, processes proceeding through the s channel are 
suppressed compared to the t-channel, and the gg partonic contributions are much smaller 
at the Tevatron than the qq ones. We take = 1 for our evaluations; however later we will 
use the known information on the Tevatron data to restrict Cd and discuss its effect to the 



lo, 
a, 




1000 



Eg 10-1" 



1000 



1,500 2000 

Az, /GeV 



2500 



3000 




1500 2000 

A„ /GeV 



3000 




1000 1500 2000 2500 

A„ /GeV 



3000 



FIG. 1: The pp 7777 cross section at y/s = 1.96 TeV as a function of Ak for = 1 and 
du = 1.1, du = 1.5 and dy = 1.9. The individual contributions from the gg and qq subprocesses 
are grouped and shown for each channel and number of unparticles exchanged. 
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FIG. 2: The pp 7777 cross section at y/s = 1.96 TeV as a function of du for Cd = 1 and Ak = 1 
TeV and Au = 3 TeV. The individual contributions from the gg and qq subprocesses are grouped 
and shown for each channel and number of unparticles exchanged. 

maximal cross section. Looking at all the contributions, the cross section is dominated by 
the qq partonic contribution with one unparticle exchange proceeding through the t-channel, 
and could reach as much as around 1 fb. The cross section is very sensitive to the values of 
the (Ik parameter, and is maximal for smaller du. 

In Fig. [2] we show the dependence of the total cross section at the Tevatron on the du 
parameter for two fixed values of the unparticle scale Ak = 1, 3 TeV. As each unparticle 
intermediate state contributes to the amplitude a factor of , the drop in cross section 

with increasing Au is expected. Thepp — > 7777 cross section at -^i = 1.96 TeV is dominated 
again by the qq partonic contribution with one unparticle exchange proceeding through the 
t-channel, and for A^ = 1 TeV, rf^^ = 1.1, we recover the maximal 1 fb) value for the cross 
section. As a common feature of both parameter investigations, the cross section for the 
process pp 7777 is completely dominated by one set of diagrams only, the t channel gq 
with one unparticle exchange. We should also note that the diagrams with two U exchange 
in the t-channel, given in the panel (c) of Fig. [11] in Appendix, are proportional the mass of 
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the quarks and we neglect them. This is also true for the subprocess qq — * 77(7(7 • 



2. The process pp — > 7799 at Tevatron 
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FIG. 3: The pp — > 779(7 cross section at ^/s = 1.96 TeV as a function of Ky for = 1 and 
du = 1.1, dii = 1.5 and djj = 1.9. The individual contributions from the gg and qq subprocesses 
are grouped and shown for each channel and number of unparticles exchanged. 



We perform the same analysis for the process pp —>■ 'j'jgg at the Tevatron. Unlike the 
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FIG. 4: The pp — > ^^gg cross section at y/s = 1.96 TeV as a function of du for = 1 and A^^ = 1 
TeV and Ajj = 3 TeV. The individual contributions from the gg and qq subprocesses are grouped 
and shown for each channel and number of unparticles exchanged. 

four photon signal, this involves two jets which would be affected by background signals 
and whose identification requires further studies. The Feynman diagrams corresponding to 
this process are given in the Appendix. As before, the relevant parameters are du and Au 
and we use the settings as in the four-photon case. In Fig. [3] we show the dependence of 
the total cross section on energy scale of the unparticle, for three values of du (1.1, 1.5 
and 1.9); and in Fig. |4]we give the corresponding dependence on the du parameter, for two 
values of Au = 1, 3 TeV. The cross section in this case is dominated by two contributions: 
the partonic contributions from qq and gg containing the minimal number of unparticle as 
intermediate states (one in this case), and without s-channel suppression (both t-channels). 
The cross section is further enhanced here with respect to pp 7777 by the presence of 
two powers of the strong coupling constant replacing the electroweak one a^, and could 
reach about 100 fb for the Au = 1 TeV, du = 1.1 case. Fig. H] shows that the maximal value 
for the cross section decreases by two order of magnitudes when we increase Au from 1 to 
3 TeV. It is interesting to note that while gg contribution is suppressed at the Tevatron for 
hadron production, and in unparticle-mediated pp 7777 , it is very strong for the process 
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PP ~^ 1199 ) while all the other subprocesses are subdominant. 
B. The multiphoton processes at LHC 
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FIG. 5: The pp 7777 cross section at ^/s = 14 TeV as a function of Au for = 1 and du = 1.1, 
du = 1.5 and du = 1.9. The individual contributions from the gg and qq subprocesses are grouped 
and shown for each channel and number of unparticles exchanged. 



We turn our attention to the four photon production and two photon plus two gluon at 
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FIG. 6: The pp 7777 cross section at y/s = 14 TeV as a function of du for = 1 and = 1 
TeV and Au = 3 TeV. The individual contributions from the gg and qq subprocesses are grouped 
and shown for each channel and number of unparticles exchanged. 

the LHC and discuss these two channels separately. We take throughout the center-of-mass 
energy ^/s = 14 TeV. We expect our signals to be significantly enhanced here, and even 
more relevant compared to the Tevatron results if the unparticle scale An turns out to be 
large. 

1. The process pp — > 7777 at LHC 

At the LHC, pp — *• 7777 can proceed through gg and qq at the partonic level. At 
^/s = 14 TeV the gluonic components of the proton dominate, so one would expect the 
cross section to be dominated by the gg contribution. This is true only for a very small 
part of the parameter space, < 1.2 TeV, and the result is dependent of the value of du- 
(For instance, for du = 1.5, the qq contribution dominates over the whole parameter space, 
whereas for du = 1.1 and 1.9 there is small region of the parameter space where the gg 
contribution is slightly larger.) The reason is that the qq partonic cross section can proceed 
through a single unparticle, while the gg requires at least two. The energy scale suppression 



13 



associated with each unparticle overwhelms the advantage from the partonic distribution 
function of gluons versus quarks in the proton, and the correlation between the number of 
unparticles in a graph and the size of its contribution to the cross section dominates the 
cross section. Here as before there is also an indication of s-channel suppression, thus the 
smallest contribution will be given by the process with three unparticles going through the 
s channel, with (in that instance) gg contribution dominating the qq one. 

In Fig. [6] we recover the expected Au suppression (three orders of magnitude going from 
Ak = 1 TeV to Ak = 3 TeV). We also note that the total relative contribution coming 
from different channels is changing: at Au = 1 TeV mainly gg s and t channels with 
two unparticle exchanges dominate, and the qq t-channel with one unparticle exchange 
contributes significantly, while at Au = 3 TeV the cross section is dominated by the one- 
unparticle qq t-channel. This is indeed expected since larger values of Au suppress further 
the channels with many unparticle exchanges. The cross section in pp —* '-f'~f'~f'~f can reach 
about 10 fb, for Au = 1 TeV and du = l-l- (Again to be expected, as the suppression 
associated with each unparticle is inversely proportional to A^".) 



2. The process pp "fygg at LHC 

Finally, in Figs. [7] and [H] we show the dependence of the pp Tigg cross section with the 
energy scale Au and dy parameter, respectively, at the LHC. The cross sections obtained in 
this case are the largest shown in this analysis, and can reach a spectacular 20 — 25 pb for 
the most advantageous parameter point Au = 1 TeV and du = 1.1. The cross sections are 
dominated by the gg t-channel partonic contribution with only one unparticle in the process, 
although perhaps slightly less so for du = 1.9. Increasing Au from 1 to 3 TeV decreases 
the cross section again by two orders of magnitude. It also slightly re-orders the relative 
contributions coming from different channels, but not significantly. 



C. The Tevatron bound on 



So far we set the three-point unparticle coup 
parameter is practically unconstrained. In Ref. 



ing constant Cd as unity. However, this 
9| an upper bound is obtained from the 



contribution of the three point interactions using the multiphoton events in the Tevatron 
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FIG. 7: The pp — > 77(75 cross section at ^/s = 14 TeV as a function of A^^ for Cd = 1 and d^^ = 1.1, 
du = 1.5 and du = 1.9. The individual contributions from the gg and qq subprocesses are grouped 
and shown for each channel and number of unparticles exchanged. 

data with integrated luminosity C = 0.83/6"^. Based on our results in previous subsections, 
we recalculate here the upper bound on Cd by including all contributions to multiphoton 
signal. The Standard Model predicts negligible background events, and no signal event is 
observed above the Standard Model prediction. Then at 95% CL the bound is 



C X app^i^ < 3.04 
15 



(14) 
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FIG. 8: The pp Jjgg cross section at y/s = 14 TeV as a function of du for = 1 and A^^ = 1 
TeV and Au = 3 TeV. The individual contributions from the gg and qq subprocesses are grouped 
and shown for each channel and number of unparticles exchanged. 

where the right-hand side is the number of quahfied events in the data. Here o"pp_»4-y = 
+ ^pf^i'y- Writing the total cross section as a sum of terms grouped in powers of 
Crf, we replace the cross section app^^^ in the left-hand side of the above expression by an 
quadratic equation 

The coefficients of various Cd terms are the corresponding cross sections calculated with 
Cd = 1 in the preceeding subsections. Note that we cannot extract the An factors from 
the terms in the left-hand side of Eq. (fT5l) . since each scales differently with A^. So, the 
upper bound on Cd can now be computed as the positive root of the above equation from 
our numerical study. The results are shown in Fig. O We separately include the results 
for the three-point interactions for comparison. Note that our cross section values from the 
three-point interactions are about a factor of 2 greater than the ones given in [9|], which 
may be due to different set of cuts used in each case as well as numerical precision in the 
computation of the function Fy, where we had to keep the number of iterations rather low. 
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It is seen that, in comparison to the bounds from the three-point interactions only, the upper 
bounds on Cd are significantly reduced when all contributions are included. The deviation 
is much bigger for smaller values of dy. The order of reduction ranges from an order of 
magnitude to more than three orders of magnitude, depending on the values of du and h^y. 
We obtain much stronger bounds for smaller djj and Ku parameters region. For example, 
the upper bound on Cd becomes around 4 at c/^ = 1.1 and Ku = 1 TeV while it can reach 
8 X 10^ if only the three-point unparticle interaction cross section is kept. Thus, our earlier 
choice of Crf = 1 turns out to be almost the maximal value at {du^Ku) = (1.1, ITeV) and 
should be considered conservative for any set of bigger {du.Ku) values. 

Now, what would be the effect of relaxing Cd = ^ and using maximal Cd on the processes 
ppij)) — ^ 7777 and pp{j)) 1199 at both colliders? First of all, this would not change 
the total cross section significantly, as long as the cross section is dominated by the single 
unparticle exchange in either s or t channel. Only if the interference terms between the 
three-point interactions and the rest of the interactions to the cross section are sizable, will 
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FIG. 9: The maximal Cd bomids from the Tevatron as a function of dy for various An on the left 
panel and as a function of Ku for various du on the right panel. The bounds are shown separately 
when only the three-point interactions are considered for comparison. AO ^ = 1 and = V2^/e 
are used. 
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there be such enhancement, which is hnear in C™"^. Of course, if we focus on the the case 
of individual contributions enhancement, the three-point unparticle interaction part will 
receive the largest enhancement since it scales as C*™"^ while in the rest there will be no 
change. A more precise analysis would first fit the curves in Fig. M to get an approximate 
relation for C^""^ as function of du when Ay is fixed. Then a search for maximal value of the 
cross section would use these functions to also obtain the relative contributions. We find it 
unnatural to go very large values, as the bounds on it most likely will become stronger 
once some data from LHC is available. As we did not include this in our evaluations, our 
total cross section values obtained are more conservative. 

IV. SUMMARY AND CONCLUSION 

Since the original suggestion of Georgi and the formulation of unparticle physics, 
viewed as a conformal theory with exact scale invariance coupled to the standard model at 
high energies, numerous phenomenological explorations of this idea have been developed. 

In the present work, we have presented a contribution to the study of implications of self- 
interactions of the low energy conformal sector. We concentrated on implications of such 
couplings on the multiphoton signal at the hadronic colliders, the Tevatron and the LHC. 
The processes ppij)) 7777 and 77(7(7 involve exchanges of one, two and three unparticles 
as intermediate states. The diagrams contributing to such processes involve unparticle ex- 
changes, as well as unparticle and vector boson exchanges. The situation is unusual, and 
evaluation of these diagrams does not resemble expressions obtained form ordinary pertur- 
bation theory, as graphs with one, two or three internal unparticle exchanges contribute to 
the same subprocess. As every unparticle vertex contains a factor of — while every SM 
particle couples to gluons (photons) with interaction strength as{cew), the evaluation of the 
cross section includes a sum over terms with different powers of coupling constants and of 
unparticle energy scale Au. Although this process is suppressed (in fact negligible) in the 
SM and can only proceed here though (at least one) unparticle exchange, the total cross 
section is dominated by the diagrams with the minimal number of internal unparticles. We 
have calculated every contribution to the cross section separately, discriminating between 
s and t channels, number of internal unparticles, and partonic components. We confirm 
earlier results for the contribution of the three unparticle exchange and correlation function. 
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but also show that this contribution is sub dominant. 

At the Tevatron and the LHC, multiphoton signals coming from unparticle exchanges 
give rise to significant cross sections, the most spectacular of which could be as large as 
20 — 25 pb, in pp — > '-fygg at LHC, for the most advantageous parameter point, the min- 
imum Ak and du investigated (A^^ = 1 TeV and du = 1.1). We consider first the case in 
which Cd, the coupling from the three-unparticle correlation function is set to one. We then 
restrict the Cd parameter based on the Tevatron data. This was done previously, assuming 
the multiphoton process is dominated entirely by the three-unparticle exchange, where ex- 
pressions for Cd can be obtained independently of the unparticle scale An. In a complete 
evaluation, terms with different numbers of unparticles enter the summation, and the re- 
striction on Cd becomes a restriction on the quadratic roots of a A^^-dependent equation, 
rather than a simple inequality. For completness, we present maximal bounds on Cd for 
both three-unparticle exchang well as for the total cross section. 

In concluding, the main point of our paper is to show that multiphoton signals at hadron 
collider are significant if one includes the coupling of SM particles with one or several un- 
particles, and if observed, these events would be a strong indication of unparticle physics. 
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VI. APPENDIX 

In this section we give all the relevant Feynman diagrams for the processes studied. We 
include only diagrams in which at least one scalar unparticle Us contributes. The Standard 
Model diagrams are not included. 
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A. Feynman Diagrams for the process pp{p) — > 7777 

We start first with the pp{p) — > 7777 and group the figures into two: one for the gg 
partonic contribution, and the other for qq partonic contributions. Following the numerical 
estimates, we group further together graphs with the same number of unparticles in the 
intermediate states, and the processes going through the s and t ~ u channels together. 
They are all given in Figs. [TOl [TT|, and [121 In all the diagrams, the wave and the spiral lines 
represent photon and gluon fields, respectively. 




12345 



FIG. 10: The Feynman diagrams contributing to the subprocess gg 7777 with the three-point 
unparticle vertex (panel a), the s-channel (panel b) and the t- and u-channels (panel c). 
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FIG. 11: The Feynman diagrams contributing to the subprocess qq — > 7777 with the three-point 
unparticle vertex (panel a) , the s-channel (panel b) and in the t- and u-channels with two unparticle 
Us exchanged (panel c). The quark qij represent all five light quarks and i = j. 
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FIG. 12: The Feynman diagrams contributing to the subprocess qq — > 7777 in the t-channel with 
a single unparticle Ug exchanged. The corresponding u-channel diagrams are also included. 
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B. Feynman Diagrams for the process pp{p) — > 7755 

Next we present the Feynman graphs for the process pp{p) — > 775'5'- The organization 
of these is the same as for pp{p) 7777: we classify contributions from different partonic 
components first, then further group them according to the number of unparticles and the 
channel through which they proceed. They are all given in Figs. [131 El [13 and [T6l 




12 3 4 



id) 

FIG. 13: The Feynman diagrams contributing to the subprocess gg 7759 with the three-point 
unparticle vertex in the s-channel (panel a), the three-point unparticle vertex in the t-channel 
(panel b), the s-channel with single unparticle exchange (panel c) and the s-channel with two 
unparticle exchange (panel d). 
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FIG. 14: The Feynman diagrams contributing to the subprocess gg — > 'j'jgg in the t-channel with 
a single unparticle Ug exchanged (panel a) and in the t-channel with two unparticles Us exchanged 
(panel b). The corresponding u-channel diagrams are also included. 
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(c) id) 
FIG. 15: The Feynman diagrams contributing to the subprocess qq — > 77(7(7 with the three-point 
unparticle vertex (panel a), the s-channel with a single unparticle Ug exchanged (panel b) and 
the t-channel with two unparticles Ug exchanged (panel c) and in the t- and u-channels with two 
unparticle Us exchanged (panel d). The quark qij represent all five light quarks and i = j. 




FIG. 16: The Feynman diagrams contributing to the subprocess gg Jjgg in the t-channel with 
a single unparticle Ug exchanged. The corresponding u-channel diagrams are included. 
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